Critical correlations in an ultracold Bose gas revealed by means of a temporal 

Talbot-Lau interferometer 
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We study experimentally the critical correlation in an ultra-cold Bose gas with a temporal Talbot- 
Lau (TL) interferometer. Near the critical temperature, we observe a bi-modal density distribution 
in an ultra-cold Bose gas after the application of the TL interferometer. The measured fraction of 
the narrower peak in the density distribution displays a clear peak within the critical regime. The 
peak position agrees with the critical temperature calculated with the finite-size and interaction 
corrections. The critical exponents are extracted from the peak and they agree with the critical 
exponents for the correlation length. 
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Near the second-order phase transition [lHij] from a 
normal fluid to a superfluid characterized by a complex 
order parameter, the diverging spatial correlation length 
around the critical temperature is the driving force be- 
hind various critical phenomena. There have been enor- 
mous experimental efforts to study the critical phase 
transition to superfluid in liquid helium One of 

the endeavors was even carried out in space to get rid 
of the deviation caused by gravity Q. Many critical 
phenomena, which had long eluded direct experimental 
study, have now been studied experimentally with ultra- 
cold Bose 0-11 1 and Fermi gases [Hf- In particular, the 
second-order phase transition from a normal state to a su- 
perfluid state has been studied in a landmark experiment 
[l3| , where the critical behavior in the correlation length 
was revealed by detecting the interference of two released 
atomic clouds. This experiment demonstrates that the 
spatial correlation in a system can enhance the interfer- 
ence effect, which in turn can be used to detect spatial 
correlation. Recently, there are also intensive theoretical 
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studies on the critical behavior with cold atoms 



In this Letter we study experimentally the spatial cor- 
relation in an ultra-cold Bose gas with a temporal Talbot- 
Lau (TL) interferometer. In such an interferometer, the 
interference effect between the correlated atoms reported 
in Ref. [l3| is greatly enhanced because all atoms in the 
Bose gas are involved. With this TL interferometer, 
we explore how the spatial correlation in an ultra-cold 
atomic gas varies as it undergoes the phase transition 
from a thermal gas to a Bose-Einstein condensate (BEC). 
At the critical temperature, we observe a bi-modal den- 
sity distribution after the application of the TL interfer- 
ometer. When the fraction of the narrower density dis- 
tribution is measured, we find a clear peak in the critical 



regime and the peak position agrees well with the critical 
temperature calculated for this Bose gas by including the 
finite-size and interaction effects. We have also extracted 
the critical exponents from the measured peak and they 
are very close to the theoretical critical exponents for cor- 
relation length, suggesting that the fraction measured in 
our experiment is proportional to the correlation length. 

The temporal TL interferometer for ultra- cold 
atoms. — The TL interferometer has been widely used 
to reveal the interference property for different systems 
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20]. Here we use a temporal TL interferometer [21 
to experimentally study the phase transition to a BEC. 
Our TL interferometer (Fig. [IJ consists of two optical 
lattice pulses separated by a time interval ry equal to the 
odd times of half a Talbot time T T = m\ 2 /2h [H [22 1 
with A being the laser wavelength, m the atomic mass 
and h the Planck's constant. In our experiment, the 
wavelength A = 852 nm, the interval ry = 3Tt/2, the 
pulse width To = 3/zs, and the lattice depth Uq = 80Er 
with En being the recoil energy of an atom absorbing 
one lattice photon. The lattice depth is calibrated 
experimentally by the Kapitza-Dirac scattering. The 
ultra-cold atoms of 87 Rb were prepared in a magnetic 
trap with axial frequency 20 Hz and radial frequency 
220 Hz [23j | . The TL interferometer was applied onto 
the ultra-cold atomic gas along the axial direction. 
The time-of-flight (TOF) images are obtained with the 
standard absorption imaging method after switching off 
the magnetic field and 30 ms free expansion. 



During the first lattice pulse, the atomic gas is divided 
into hundreds of disk-shaped atomic gases that are phys- 
ically separated due to the strength of the pulse. Af- 
ter the pulse, these disk-shaped atomic gases start to 
merge and interfere for a given time. The second lat- 
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FIG. 1: (color online) Effects on a pure condensate and a 
thermal gas by the Talbot-Lau interferometer, where two op- 
tical lattice pulses are separated by an interval r/ equal to 
the odd times of half a Talbot time. The two-dimensional 
experimental TOF images are for (al) a pure condensate and 
(a2) a thermal gas before the first lattice pulse; the momen- 
tum distribution of (bl) a BEC and (b2) a thermal gas after 
the first pulse; the momentum distribution of (cl) a BEC and 
(c2) a thermal gas after the second pulse. 



tice pulse again divides it into hundreds of disk-shaped 
atomic gases, leading to a second interference. For an 
initial atomic gas without correlation, the interference is 
in fact for each individual atom itself. The overall inter- 
ference pattern is a simple addition of the interference 
intensities for all the individual atoms in the gas. Such 
a simple addition does not apply for an atomic gas with 
spatial correlation larger than the spatial period of the 
pulsed lattice. Due to the correlation, one should instead 
sum up first the interference wave functions for the cor- 
related atoms within the correlation length, then square 
the amplitude of the summed wave function to obtain an 
enhanced interference. This enhancement suggests that 
the TL interferometer would be powerful for detecting 
the spatial correlation among atoms. 

Our TL interferometer is effectively an interferometer 
for atoms around zero momentum and with spatial corre- 
lation larger than the spatial period of the pulsed optical 
lattices. To clearly see this, we turn to the momentum 
space, where the experimental observation is made, and 
focus on atoms with momenta very close to zero. For 
these atoms, after the first optical lattice pulse, signifi- 
cant proportion of atoms will be transferred to the mo- 
mentum around ±2nhk (k = 2tt / X, n = 1, 2, 3, • • • ) with 
almost no atoms left at zero momentum; after the second 
optical lattice pulse, the atoms with momentum around 
±2nhk will be brought back to the momentum around 
Ohk [24|. This effect is verified and illustrated experi- 
mentally with a pure condensate in the left column of 

Fig.m 

With this in mind, it is straightforward to see how our 
interferometer will affect a thermal atomic cloud without 
correlation: a very small fraction of atoms that are in the 
state of momentum zero are similarly transferred away 
and brought back to momentum zero by our interferom- 
eter; for majority of the atoms with non-zero momenta, 



they will not be brought back to their original momen- 
tum states. The end result is an atomic cloud with much 
wider distribution in momentum space (see the right col- 
umn of Fig. []}. For thermal atomic cloud, because of the 
lack of correlation between different momentum states, it 
is expected that there is a trivial mapping before and af- 
ter the TL interferometer. Both experiment and theoret- 
ical simulation show that the momentum distribution is 
still a Gaussian distribution after the TL interferometer. 

Far below the critical temperature which comprises 
both condensate and thermal atoms, after the applica- 
tion of the TL interferometer, we expect a bi-modal den- 
sity distribution because the condensate will still lead to 
a narrower central peak as experimentally shown in Fig. 
Hb5). 

When the system is close to the critical temperature, 
there is strong spatial correlation between atoms, which 
will enhance the interference as reported in Ref.[13[. As 
our TL interferometer is essentially an interferometer 
for atoms around momentum zero, this enhancement is 
expected to increase the population around momentum 
zero, causing the momentum distribution around zero to 
deviate from the Gaussian distribution after the applica- 
tion of the TL interferometer. Such a deviation from the 
Gaussian should be bigger for stronger correlation. In a 
sense, the existence of an order parameter for the corre- 
lated atoms within the critical regime makes these atoms 
have similar behavior to a condensate after the applica- 
tion of the TL interferometer. This physical picture will 
lead to a bi-modal density distribution after the applica- 
tion of the TL interferometer, similarly to the system far 
below the critical temperature. 

Within the critical regime marked by the Ginzburg 
temperature Tq. i.e. \T— T c \ < Tq, the correlation length 
diverges near the critical temperature T c as [H-Q 
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where the critical exponents v = vt = 0.67. This is for 
the infinite system in the thermodynamic limit. For the 
finite system in our experiment, the correlation length no 
longer diverges. Nevertheless, the correlation is still the 
strongest at the critical temperature T c and the correla- 
tion length should have a peak around T c . Within the 
critical regime, this implies that the deviation from the 
Gaussian distribution after the TL interferometer is the 
largest at T c . 

Observation of the critical phase transition. — With 
the above expectations, we have surveyed the ultra-cold 
atomic gases with the TL interferometer over a wide 
range of temperatures with particular attention paid to 
the range where the critical phase transition from a ther- 
mal cloud to a BEC occurs. The results for five typical 
temperatures are shown in Fig. O The system tempera- 
ture is expressed in terms of Tj? = O.QAhLJhoN 1 ^ 3 /ks, the 
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critical temperature of the corresponding ideal Bose gas 
without the finite-size correction [25J ■ u>ho = (^x^y^z) 1 ^ 3 
is the geometric average of the harmonic trapping fre- 
quencies and N is the total particle number. 

Shown in the first row of the figure are TOF images 
before the application of the TL interferometer. The 
corresponding integrated density distributions are shown 
as dotted blue lines in the third row. There is a clear 
bi-modal distribution at temperature T/T c ° = 0.7 (Fig. 
|2](a5) ) and a minor bi-modal distribution at tempera- 
ture T/T° = 0.89 (Fig. H>4)). There are no obvi- 
ous condensate fraction at three other higher temper- 
atures. However, after the application of the TL in- 
terferometer, a clear bi-modal (i.e., non-Gaussian) dis- 
tribution emerges at T/T c ° = 0.92 (Fig. [2£b3)). In 
contrast, at a temperature well above the critical tem- 
perature (T/Tj? = 1.1), and temperatures slightly lower 
(T/T c ° = 0.89) and slightly above (T/T c ° = 0.95), there 
are no obvious bi-modal structures after the application 
of the TL interferometer. This shows that the emergence 
of the bi-modal distribution at T/Tj? = 0.92 is the re- 
sult of the diverging correlation length near the critical 
temperature. 

To quantify the observed critical behavior, we use the 
fraction of atoms in the central peak of the bi-modal 
structure after the TL pulses for different temperatures. 
For each density distribution obtained after the TL inter- 
ferometer, we fit it along the direction of the TL pulses 
with a bi-Gaussian function, and compute the fraction 
of the central narrow peak with f r = 1 — A^/At, where 
Ab is the area under the broad peak and At is the total 
area in the bi-modal structure. Note that it is an as- 
sumption that the non-Gaussian distribution of momen- 
tum observed near T c is bi-modal. This assumption is 
reasonable as indicated in Fig. [U(b3) (see Supplemental 
Material for further discussions). 

In Fig. [31 we show how this fraction f r changes with 
temperature. Near the critical temperature, a clear peak 
around T/T° = 0.92 is seen. We also notice a small dip in 
the peak, which may due to the uncertainty of the tem- 
perature calibration which is about 0.011"°. Because the 
central peak after the TL interferometer physically orig- 
inates from the spatial interference of the atoms, within 
the critical regime, it is expected that the larger the cor- 
relation length, the larger the fraction f r . Hence, this 
peak implies the critical phase transition from a thermal 
cloud to a BEC. 

As the correlation length diverges at the critical tem- 
perature, this peak position can be regarded as the crit- 
ical temperature of the system. Above this critical tem- 
perature, the fraction f r increases with decreasing tem- 
perature, very similar to what was observed with the 
method of the matter- wave interference [l3[ • Below this 
critical temperature, we see f r decreases with decreasing 
system temperature within the critical regime. 

From the peak in Fig. |31 we arc able to determine the 



critical temperature T c . We find T c /T c ° = 0.92 ± 0.01, 
indicating a negative shift from the critical temperature 
T c ° of the ideal Bose gas. The negative shift is due to the 
finite-size and interaction corrections of the critical tem- 
perature. The simultaneous consideration of the finite- 
size and interaction corrections to the critical tempera- 
ture [26, 27 1 gives a negative shift of about 0.08 agreeing 
well with our experimental result T c /T c ° = 0.92 ± 0.01. 
The correction to the critical temperature also agrees 
with previous experiments 29|-32j|. These agreements 



further confirm the observation of the critical phase tran- 
sition in this work. It is worth to point out that we have 
computed the Ginzburg temperature Tq/T® J26|. The 
result is that \T G - T c |/T c ° w 0.06, which agrees with the 
peak width in Fig. G3 

Blow the critical temperature and within the criti- 
cal regime, there is a decrease of f r . This is because 
that the fraction f r has two origins, the condensate frac- 
tion that has infinite coherent length and the correlated 
thermal fraction that obeys Eq.([2]). Within the criti- 
cal regime, the contribution from the correlated thermal 
fraction dominates, causing the fraction f r to decrease. 
In fact, we measured the condensate fraction without the 
TL interferometer. The results agree with the theoretical 
result by taking into account the finite-size and interac- 
tion corrections [26T-|2l|. At T = 0.85T C °, just outside the 
critical regime and where the fraction starts to increase, 
our experiment shows that the condensate fraction with- 
out interferometer is about 0.08 which is of the order of 
/,■ in Fig [3l We emphasize that this understanding of the 
whole system consisting of two parts, condensate and cor- 
related thermal atoms, is phenomenological in the spirit 
of two-fluid model for liquid helium. A more rigorous 
theory is needed to fully understand the strongly corre- 
lated quantum many-body state just below the critical 
temperature and how it responds to an interferometer. 

Discussion and outlook. — Within the critical regime, it 
is clear that the fraction f r is a monotonic function of the 
correlation length. It is natural to fit the data on the right 
side of the peak with a/{T/T° - T c /T°f + b and the left 
side with a'/(T c /T c - T/T^f +b' + c'(T c /T° - T/T°). 
The additional term for the left side is to account for 
the condensate fraction which appears below the crit- 
ical temperature. Rigorously speaking, the correlation 
length should display a smooth evolution within the crit- 
ical regime due to the finite-size characteristic of our ex- 
periment. Similarly to Ref. [l3j], because of the uncer- 
tainty of the temperature calibration, we do not discuss 
the finite-size effect to the correlation length in this work. 

From the data fitting, we get v' = 0.70 ± 0.08 and 
v = 0.70 ±0.11. The uncertainties correspond to a 95% 
confidence level. These two critical exponents are very 
close to v = vt = 0.67, the theoretical value of the corre- 
lation length's critical exponent for the universality class 
of three-dimensional XY model [33j, |34[. If we assume 
that f r ~ and the theoretical value of the critical 
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FIG. 2: (color online) Revelation of the critical phase transition from a thermal cloud to a BEC by the TL interferometer. 
The first row shows the density distributions of the atomic gases before the TL interferometer for five different temperatures 
T/T° = 1.10, 0.95, 0.92, 0.89, 0.70. The density distributions after the TL interferometer are in the second row. The dotted blue 
lines and dotted red lines in (cl)-(c5) give the one-dimensional density distributions after the integration of the two-dimensional 
density distributions along the vertical direction for (al)-(a5) and (bl)-(b5), respectively. A bi-modal structure emerges after 
the TL pulses in (b3) and reveals the critical phase transition occurring around T/T c ° = 0.92. In (cl)-(c5), the solid black lines 
give the Gaussian fitting to the broad peak of the density distribution after the TL interferometer. 



exponent is exact, we have a ~ 1.04 ± 0.12. This experi- 
mental result indicates that the fraction f r near the peak 
in Fig. |3] is proportional to the correlation length. 

This encouraging indication is supported by a phe- 
nomenological theory (see Supplemental Material). How- 
ever, we still do not have a rigorous theory how a strongly 
correlated system responds to an interferometer as the 
action of an interferometer can no longer be regarded 
as a small perturbation. This is different from known 
measurements on critical systems, for example, the heat 
capacity measurement on liquid helium, which can be re- 
garded as a perturbation and have negligible effects on 
the systems. 

In summary, we have studied the critical behavior of 
interacting ultra-cold Bose gases with a temporal TL in- 
terferometer. A peak in the fraction of the narrower peak 
of the density distribution was clearly observed across the 
phase transition from a thermal gas to a BEC. This ex- 
perimental study of the ultra-cold atomic gas within the 



critical regime opens the way to relevant theoretical stud- 
ies of critical dynamics, which has many open problems 
[H-Q. It also gives a new method to measure the critical 
temperature. We expect that the temporal TL interfer- 
ometer be used to study a wide range of the critical phase 
transitions, such as the quantum phase transition from 
supcrfluid to Mott insulator for cold atoms in an opti- 
cal lattice [H , ultra-cold Fermi gas [l2[ and the quantum 
magnetism with the cold atoms and molecules as a quan- 
tum simulator 
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FIG. 3: (color online) The measured fraction f r of the central 
peak as a function of temperature. Each point is an average 
over about ten experimental data within a temperature step 
5T/T® = 0.01. The two solid lines are fittings to the data near 
the peak with Eqs. (1) and (2). The colored areas represent 
the uncertainties with a 95% confidence. The vertical dash 
line shows the location of the theoretical critical temperature. 
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Supplemental Material 



Temporal Talbot-Lau interferometer for a condensate. Our temporal TL interfer- 
ometer has been described in the main text and its effects are illustrated in Fig. 1 of the 
main text. We here provide a detailed theory and explain why the momentum distribution 
of the condensate evolves as shown in the left column of Fig. 1 of the main text under the 
TL interferometer. 

Consider a condensate with very narrow momentum distribution around zero. After the 
first short optical lattice pulse, which is described as V(z) = Uocos 2 (kz), the atoms are in 



a state approximately described by [SI] 



(to) oc exp{i ^° T ° cos 2 (£;,z)} |0) 
a 



oo 

OC 

n=— oo 



i n Jn(U T /2h)\2nhk), (S.l) 



where \2nhk) (n = 0, 1, 2, ■ • • ) represents a state of momentum 2nhk and J n is the Bessel 
function of the first kind. Our experimental parameters have been chosen such that 
Jo(UoTo/2h) = 0, which means that the zero momentum is no longer occupied after the 
first pulse (Fig. l(bl)). 

After a free evolution of time 77, the second pulse which is identical to the first one is 
applied. The quantum state is again transformed and becomes 

00 00 

^(27-o + T,) oc iU [ J2 Jn-Mr /2h)e- ie -J m (U r /2h)]\2nhk), (S.2) 

n=— 00 m=— 00 

where 9 m = 2irm 2 Tf /Tt- Since r, = 3T?/2 in our experiment, we have 

00 

0(2r o + 3T T /2) oc ^ i n J n (0)\2nhk) = |0). (S.3) 

n=— 00 

This means that the second pulse has completely brought the system back to the state of 
zero momentum (Fig. l(cl)). 

Note that in getting Eq. ( 1S.2j) . there is an implicit assumption that during the second 
optical lattice pulse, there is full spatial overlapping between different momentum state 
\2nhk). For our experimental parameters for a pure condensate, this condition holds. Also 
note that in the above derivation, we have ignored the non-homogeneous distribution of the 
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condensate in a real experiment. It is clear from the experimental result in Fig. 1 that this 
non- homogeneity is not important. This is also verified by our numerical simulation. 

The effect of our TL interferometer is observed experimentally with a pure 87 Rb BEC as 
shown in the left component of Fig. 1, where the experimental time-of-flight (TOF) images 
are given. The TOF images are obtained with the standard absorption imaging method after 
switching off the magnetic field and 30 ms free expansion. A cigar-shaped pure condensate 
of 10 5 atoms was created with axial frequency 20 Hz and radial frequency 220 Hz. Its 
momentum distribution is shown in Fig. l(al). When only the first pulse is applied along 
the axial direction, as shown in Fig. l(bl), most of the atoms (about 60%) are transferred 
to the momentum around ±2hk. About 36% of the atoms are populated around ±4hk with 
the population around momentum zero almost completely depleted. After a time interval of 
37V/2, the second pulse with the exact same parameters is applied, as shown in Fig. l(cl). 

For thermal atomic cloud far away from the critical temperature, because it is an 
incoherent superposition of different momentum states, after our TL interferometer, we 
observe a broadened momentum distribution shown on the right of Fig. 1. Our numerical 
simulation also gives this behavior. 

Temporal Talbot-Lau interferometer for the system near the critical tempera- 
ture. In our experiment, the TL interferometer is applied along the axial z direction of 
a cigar-shaped atomic cloud, and we can treat it as an effectively one dimensional system. 
Near the critical temperature, we can introduce an order parameter ip (z) = (z)) to 
describe the system based on the widely used phenomenological theory for phase transition. 

The evolution of the order parameter ip (z, t) can be studied with the following propagator 
method 



Here K (z, t; Zi,to) is the propagator, describing how the system evolves. For our experiment, 
one may regard it as the evolution caused by the application of the TL interferometer. The 
density distribution is then 



n (z, t) = ip* (z, t) ip (z,t) = / K* (z,t;z 1 ,t )K (z,t;z 2 ,t )ip* (z 1 ,t )ip(z2,to)dz 1 dz 2 . (S.5) 



Within the critical regime, the order parameter has significant amplitude and phase 
fluctuations. We can write the order parameter as ip(z,t ) = ipo (z,t ) 9 (z, t ). Here 




(S.4) 
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ipo(z,t ) = d(W(z,to)V(z,to)) is a smooth function with uniform phase while 9(z,t ) 
describes the amplitude and phase fluctuations. With this, Eq. (1S.5j) becomes 



n (z, t) = I K* (z, t; z l} t ) K (z, t; z 2 , t ) ipo (zi,t ) 9* (z 1 ,t ) ip (z 2 , t ) 9 (z 2 , t ) dz l dz 2 . 

(S.6) 

U)) = G(z 1 ,z 2 ), where 



Z-2 



S2 



S3j | . After the ensemble 



We now take the ensemble averaging, which yields {9* (zi,t ) 9 
G (zi, z 2 ) is the correlation function between two points Z\ and z 2 
averaging, we obtain the population of the atoms with the second-order correlation, 

n (z, t) = f K* (z, t; z u t ) K (z, t; z 2 , t ) ip* Q {z u t ) (9* (z^to) 9 (z 2 , t ))ip {z 2 , t ) dz 1 dz 2 

= J K* (z, t] Zi,t Q ) K (z, t; z 2 , t ) ipQ {zi,t ) G (zi, z 2 ) tp (z 2 , t Q ) dz x dz 2 . (S.7) 

Although the above expression is rather general and applies to almost all systems, useful 
information can still extracted. When the system is far above the critical temperature, we 
have G = and n = as expected. When the system is far below the critical temperature, 
where it is dominated by the condensate fraction, the correlation function is G(zi,z 2 ) = 
1. Hence, hq relates directly to the condensate. In general, it is clear from the above 
equation, the population of the correlated atoms is a monotonic functional of the second- 
order correlation function G (zi, z 2 ): when G is strong, the population n is larger. This fact 



can qualitatively explain the experimental observations in Ref. [S4j and our experiment. 
However, the exact relation is complicated by the existence of other functions, such as the 
propagator K, in Eq. (1S.7j) . As it is difficult to obtain the exact form of K for this strongly 
correlated system, the quantitative explanation of these two experiments is elusive. 
To make progress, we write the whole density distribution n a u (z, t) as 

n a u (z, t) = n (z, t) + n incoherent (z, t) . (S.8) 

For these correlated atoms, the width of momentum distribution of n (z, t) is assumed to 
have the same as that of a pure condensate with the same form of the order parameter. 
The incoherent component, however, has a wider density distribution because of the lack 
of spatial interference after the second pulsed optical lattice. Hence we expect a bi-modal 
density distribution within the critical regime as observed in our experiment. In the main 
text, f r = J no (z,t) dzj J n a u (z,t) dz is the measured fraction of the correlated atoms. 

For the TL interferometer, the above expression can be further simplified with a discrete 
method. After the first optical lattice pulse, the atomic cloud is broken up into hundreds of 



subsystems. If the atoms in the cloud are strongly correlated, each of these subsystems will 
contribute an order parameter with a complex amplitude ipj (j denotes the j'th subsystem) 
to the total amplitude. After the second pulsed optical lattice, we have at the end of the 
TL interferometer action 

n (z, *) = 5 ^ ^ *) Gi ^i & *) ' ( S - 9 ) 

ij 

where Gij = G(\zi — Zj\) is the correlation function between two points Zi and Zj. 

During the interference process, each subsystem ipj(z,t) will expand and become almost 
completely overlap with each other. At the same time, within the critical regime, the 
correlation length £ in G(\zi — Zj\) is much larger than the spatial period of the pulsed 
optical lattice while much smaller than the length of the system (although we can not 



measure directly the value of the correlation length, the work in Ref. [S4j ] with similar 
parameters shows that the correlation length in the critical regime satisfies this condition). 
With all these consideration, we have 

n (z, t) oc X ( z i t) ipj 0, t) X G ij X ^ ^' ^ ^' *) / / dzidz&GQzi - z 2 \) . 

ij ij ij 

(S.10) 

Another reason that the above approximation is good lies in that the atoms concentrate in 
the central region of the magnetic trap, thus the homogeneous approximation can be used. 
In the critical regime above the critical temperature, the correlation function is 

G{z) oc eXp( 7 N/0 (S.ll) 
\z\ 

temperature, the correlation function is dominated 



In the critical regime below the critica. 
by the transverse part and we have S5] 



G{z) oc pr. (S.12) 

\z\ 

Nevertheless, after a simple integral transform, we have for both cases 

n oc £ . (S.13) 

Therefore, we have the fraction of the correlated atoms f r oc £ in the critical regime. This 
simple relation shows why the TL interferometer has the ability to reveal the critical behav- 
ior. 



4 



The above analysis relies on two assumptions. (1) The first pulsed optical lattice will 
not change the correlation function Eq. (IS. Ill) or Eq. (IS.12j) when the correlation length is 
larger than the spatial period of the pulsed optical lattice. This is a valid assumption for 
the experimental pulse width To = 3/is, which is much shorter than the thermal equilibrium 
time. (2) There is no thermal equilibrium between two pulsed optical lattices in the TL 
interferometer. In our experiment, the interval between two pulsed optical lattices is Tf = 
37V/2, which is much smaller than the thermal equilibrium time of the order of several 
milliseconds. 

With the above analysis on the correlated atoms, we see that the TL interferometry has 
the effect of separating the correlated atoms and incoherent atoms. Because the correlated 
atoms will be brought back to the state near zero momentum after the TL interferometry, 
the density distribution of the correlated atoms leads to the central peak in the whole density 
distribution. This leads to the bi-modal structure in the density distribution after the TL 
interferometer, if there is significant fraction of correlated atoms. 

For T/T c ° = 0.70, without the application of the TL interferometer, we see a clear bi- 
modal structure which shows the density distribution of the condensate and thermal atoms, 
respectively. After the application of the TL interferometer, there is no essential change in 
the bi-modal structure, because the condensate will still appear in the central peak in the 
bi-modal density distribution, while the thermal atoms give only the wider background. At 
the critical temperature T/T c ° = 0.92, however, the bi-modal structure appears after the 
application of the TL interferometer. The only interpretation is that the correlated atoms 
will appear in the central peak although there is no condensate in this situation. 
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